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Abstract 

, of state w < —1, we investigate the cosmological evolution of such a phantom field in a specific 

^ ' potential, exponential potential in this paper. The phase plane analysis show that the there is a 

late time attractor solution in this model, which address the similar issues as that of fine tuning 

(N ; 

. problems in conventional quintessence models. The equation of state w is determined by the 



attractor which is dependent on the A parameter in the potential. We also show that this model 



> 

Q ' is stable for our present observable Universe. 
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1. Introduction 



Recent observation shows that our universe is made up of roughly two third of dark 
energy that has negative pressure and can drive the accelerating expansion of the universe^]. 
Present data from the observation allows the equation of state in the range —1.62 < w < 

" n 

— 0.74[2]. However, the equation of state of conventional quintessence models j3|that based 
on a scalar field and positive kinetic energy can not evolve to the the regime of w < — 1. 



Some authors I 



eld ana positive Kinetic energy can not ev 



investigated a phantom field model 



which has negative kinetic energy and can realize the < — 1 in its evolution. Although 
the introduction of a phantom field causes many theoretical problems such as the violation 
of some widely accepted energy condition and the rapid vacuum decay Q], it is still very 
interesting in the sense that it can fit current observations well. Comparing with other 
approaches to realizing the w < —1, such as the modification of Friedman equation, it seems 
more economical. 

In this paper, provisionally leaving aside the theoretical puzzles about phantom field, 
we study the detailed evolution of the phantom field and the attractor property of its so- 
lution via a specific model, by using the qualitative approach to the dynamical system, 
phase plane analysis, we prove the existence of a late time attractor solution, at which the 
phantom become dominant and the equation of state is freezed with only small oscillation. 
Exponential potentials attract much attention because they can be derived from the effective 
interaction in string theory, higher dimensional gravity and the non-perturbative effects such 
as gaugino condensation [l^ and their roles in cosmological context have also been widely 
investigated 13]. We therefore consider the phantom evolution in exponential potential in 
the following. 

2. Phantom in Exponential Potential 

In the spatially flat Robertson- Walker metric, 

ds"^ = dt^ - (t) {dx^ + dy"^ + dz^) (0.1) 
The Lagrangian density for the spatially homogeneous phantom field is 

L^ = -\g'''{d,4>)m)~V{(t>) (0.2) 
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when consider the presence of baryotropic fluid, the action for the model is 

S= I d^x^i-^Rs -p, + L^) (0.3) 

where g is the determinant of the metric tensor g^i,, — SttG, Rg is the Ricci scalar, and pj 
is the density of fluid with a baryotropic equation of state P7 = (7 — 1)P7, where < 7 < 2 
is a constant that relates to the equation of state by to = 7 — 1. By varying the action, one 
can obtain the Einstein equations and the equations of motion for the scalar fleld as 

^2 



where 



P7 = -3H(pj+p^) 

^ + 3H<^ + AkF (0) = 
^2 

= y[P7 + /5*] 



(0.4) 
(0.5) 
(0.6) 
(0.7) 

(0.8) 



= - V{4>) (0.9) 

are the energy density and pressure of the $ field respectively , and H is Hubble parameter. 
The potential V{(j)) is exponentially dependent on (f) as V{(f)) — Vq exp(— AK0).The equation- 
of-state parameter for the phantom is 

™=^=e±^ (0.10) 

It is clear that the phantom field could realize the equation of state w < —1 when 

< 0^ < 2^(0) (0.11) 
3. Attractor Solution of Phantom Model 



In this section, we investigate the global structure of the dynamical system via phase 
plane analysis and compute the cosmological evolution by numerical analysis. Introduce the 
following variables 



X = -^R (0.12) 
N = log a 

the equation svstem ()U.4ll0.7|) could be reexpressed as the following autonomous system: 

^ = ^^[7(1 + y') - 2x'] - {3x + J^Xy') (0.13) 



dy _3 
dN ~ 2 



l/[7(l + x^- y^) - 2x^] - ^Xxy (0.14) 



Also, we have a constraint equation 



where 



^^ + ^ = ^ (0.15) 



^^ = ^ = y"'^" (0-16) 



The equation of state for the scalar fields could be expressed in term of the new variables 

as 

P4, x^ + y^ 

w^ = — = — (0.17) 

x^ - y^ 

It is not difficult to find the physically meaningful critical points {xc,yc) are (0,0) and 



(— y 1 + ^). To gain some insight into the property of the critical points, we write the 
variables near the critical points (xc, yc) in the form 



x = Xc + u (0.18) 

y = yc + v 
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where u, v are perturbations of the variables near the critical points. Substitute the expres- 
sion ()0.18|) into the autonomous system ()0.4ll0.7p . one can obtain the equations for the linear 
perturbations up to the first order as following: 



du 

— = Muu + Muv (0.19) 
du 

— = M21M + M22V 
dN 

The coefficients of the perturbation equations form a 2 x 2 matrix M whose eigenvalues 
determine the type and stability of the critical points. One can easily find that the eigenval- 



ues of M for (0,0) is (f (-2 + 7), ^) and those for (-^, \ l + is (-3 - ^, -87 - A 



\2\ 



Therefore, point (0,0) is a saddle point while (— y t) ^ stable node, which corre- 
sponds to an attractor solution. At this attractor solution, from Eq. ()0.16| 10.171 ). we know 
that ^0 = 1 and to^ = — 1 — which indicate the phantom domination and the possibility 
oi w < —1. 

Next, we study the system by numerical analysis. The results are in Fig.l, Fig. 2 and 
Fig. 3. The computation was done at A = 1.2 and 7 = 1.0. From the figures, we can 
find that the attractor property of the solution and the phantom field will slightly oscillate 
around the fixed point at late time and the equation of state is smaller than —1. 

4. Discussion 

We studied a specific phantom model and show that there exists a late time attractor 
solution in the evolution of the field. The late time attractor solution corresponds to the 
phantom dominate phase and the equation of state could be smaller than —1. This is a very 
interesting feature of the model we considered here. Now, we consider the stability of the 
phantom field under perturbation in such a potential. The perturbed metric in synchronous 
gauge could be expressed as 

ds^ = df - a'^{t){6ij - hij)dx'dx^ (0.20) 
Then, a Fourier mode of the phantom field 

50(t,k) = / 6(j){t,x)e-'^-''d^x (0.21) 



5 



FIG. 1: the phase graph of y to x 




FIG. 2: The evolution of Equation of state of Phantom w with respect to N 
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FIG. 3: The evolution of with respect to 



satisfies the equation of motion 



(0.22) 



where the h is the trace of hij and the prime denote the derivative with respect to (p. The 
effective mass for the perturbation is [k'^ — When the potential is chosen as the 

exponential V{(f)) = Voexp(— Ak0), we have [k"^ — A^k^Vq exp(— Ak0)]^/^. When the field 
evolves to its stable attractor solution, the cutoff wave number of the perturbation should 

be , 

k < XH\ 3 + — (0.23) 



so that the instability does not appear. Surely, it is not to say that the instability will not 
appear during the evolution before the field reaches the attractor. Since we do not know 
what will be the specific evolution track the field took in the past, we can only say that at 
present situation, which is the attractor solution epoch, the perturbation should not violate 
Eq. ()0.23p to safeguard the stability of the phantom field. Another remarkable feature of 
such an model is that the evolution of field is toward to its attractor no matter what is its 
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initial value. Since w^p = —1 — 3-, it is determined by the parameter A of the potential and 
is independent of the choice of the initial value of the field, which make the fine tuning of 
the field unnecessary. 
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